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Abstract. In a recent paper [1] Carinena et al analyzed a non-polynomial one- 
dimensional quantum potential representing an oscillator which they argued was 
intermediate between the harmonic and isotonic oscillators. In particular they proved 
that it is Schrodinger soluble, and explicitly obtained the wavefunctions and energies 
of the bound states. In this paper we show that these results can be obtained much 
more simply by noting that this potential is a supersymmctric partner potential of 
the harmonic oscillator. We then use this observation to generate an infinite set of 
potentials which can be exactly solved in a similar manner. 



1. Introduction 

In a recent paper in this journal, Carinena et al [1] investigated the solutions to the 
eigenvalue problem for the Schrodinger equation [2] 



+ 



2 .2x^-1 



(2x2 + 1)2 



^ = Eijj. 



Using some rather involved mathematics they were able to show that the eigenfunctions 
are given by 

^"^"^^^ (2^2?!) ""^'^'' where n = 0, 3, 4, 5 . . . , 
and the polynomial factors Pn{x) are related to the Hermite polynomials by 

1 n=0 

Hn{x) + 4nHn-2{x) + 4n(n - 3)iJ„_4(x) n=3,4,5. . . 
The corresponding eigenvalues are given by 

En = —3 -|- 2n where n = 0, 3, 4, 5 . . . 

In this paper we show that Carinena's potential is a supersymmctric partner potential 
of the harmonic oscillator. This allows us to rederive all the above results in a much 
simpler fashion. In addition, by considering the question of which other potentials 



Pn(x) 



Factorization solution of Carinena 's quantum nonlinear oscillator. 



2 



are supersymmetric partners of the harmonic oscillator, we construct an infinite set of 
exactly soluble potentials, along with their eigenfunctions and eigenvalues. 

The rest of the paper is organized as follows. In section 2 we give a brief summary of 
the ideas of the factorization approach to the Schrodinger equation and supersymmetric 
partner potentials. In section 3 we prove that Carinena's potential is a partner potential 
of the harmonic oscillator, and use this to rederive the formulae for its eigenfunctions and 
eigenvalues. In section 4 we show how to find other partner potentials of the harmonic 
oscillator, and derive an infinite set of exactly soluble potentials. 



2. The factorization approach 



In this section we provide a self-contained introduction to the factorization approach, 
and the idea of supersymmetric partner potentials. More details can be found in 
references [3] and [1]. 

We start with the one- dimensional single-particle Schrodinger equation. 



HMx) 



+ VAx) 



The idea is to factorize the Hamiltonian operator, ifi, which is a second-order differential 
operator, into a product of two first-order differential operators. 



Hi 



where 



d 



A = ^ + W{x) 
dx 



At = — ^ + W{x) 
dx 



Upon direct substitution we see that 

V^{x) = W\x) - W'{x). 

We now define the operator H2 = AA'^ by reversing the order of A and A^". Simple 
algebra shows that H2 is a Hamiltonian corresponding to a new potential V2{x), 

rf2 



Ho = AA^ 



+ V2{x), V2{x) = W{xY + W'{x). 



The potentials Vi{x) and V2{x) are known as supersymmetric partner potentials. 

The key result we will be using in later sections is that the eigenvalues and 
eigenfunctions of Hi and H2 are related. Suppose that ^l^^ is an eigenfunction of 
Hi with energy eigenvalue Then Atp!^^ is an eigenfunction of H2 with energy 

eigenvalue E^^ since 



Ho 



Ai). 



(1) 



AA^A^l), 



(1) 



A 



Hit 



(1) 



A 



E^^^^^A = E^^ [A^^ 



,(1) 



The only way that A-?/;^^^ could fail to be an eigenfunction of H2 is if it equals zero, 
in which case W{x) = ——Inip^^x), which could be true for at most one value of n. 
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Similarly suppose that ipn"^ is an eigenfunction of H2 with eigenvalue E^^ . Then A'^ip^^^ 
is an eigenfunction of Hi with energy eigenvalue E^^ since 

The only way that A^ip^^^ could fail to be an eigenfunction of Hi is if it equals zero, in 

which case W{x) = — In iIjI^\x), which again could be true for at most one value of n. 

ax 

It follows that, with the possible exception of one level. Hi and H2 have the same 
energy spectra. The eigenfunctions of H2 can be found from those of Hi by applying 
the operator A; the eigenfunctions of Hi can be obtained from those of H2 by applying 
the operator A^. If Hi has one more energy level than H2, this energy eigenfunction 
will be the solution of Aip^'^ — 0, and have energy eigenvalue zero; if H2 has one more 
energy level than Hi, this energy eigenfunction will be a solution of A^i/jq'^ = 0, and 
have energy eigenvalue zero. The upshot of all this is that, if we know how to exactly 
solve one of Hi or H2, we can immediately derive the exact solution of the other. 

Finally we note that wc can easily normalize the newly obtained eigenfunctions. If 
is a normalized eigenfunction of Hi with energy E^^\ the corresponding normalized 
eigenfunction of H2 is 

where C„ is a constant to be determined. If we take the scalar product of this 
eigenfunction with itself we get 

This leads to the result 



and we similarly find that 



-1/2 



(1) 



-1/2 



(1) 



',(2) 





-1/2 







-4- + wi^) 

ax 



(2) 



3. Solution of the Carinena potential 

We can solve the Carinena potential by noticing that it is a partner potential of the 
harmonic oscillator. To see this we simply take 

from which trivial calculus gives us 

Vi{x) = W{xf - W'{x) =x^ + 8 ^ "-^j + 3 



(2x + 1)^ 



V2{x) = W{xf + W\x) = + 5. 
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If follows that 



H^ = A^A = Hc + 3 = + 3 



H2 = AA^ = Hh + 5 



(2) 



Ej! + 5, 



where 



c 



. (23:^ - 1) 

Cix2 ' ' (2x2 + 1)2 

is the Carinena Hamiltonian, and 



H 



H 



is the harmonic oscillator Hamiltonian. 

Now the eigenvalues and unnormalized eigenf unctions of the harmonic oscillator [5] 
are given by 



H 



2k + 1 



and 



X] 



H,{x)e-^"'\ 



where k = 0,1,2... and the Hermite polynomials Hk{x) are given by the Rodrigues 
formula 



HJx) 



■1)V 



dx^ 



If follows trivially that the eigenvalues and unnormalized eigenf unctions of H2 are 



E 



(2) 



2k + Q 



and 



^ljf\x)=H,{x)e-^'l\ 



The corresponding eigenvalues of Hi are therefore E)^ 
eigenfunctions 



(1) 



d 4x 

\- X -\ 

dx {2x + 1) 



2k + 6, with unnormalized 



HJx)e-^"/^. 



We can simplify this expression by repeatedly using the Hermite polynomial identities 

H',{x) = 2kHk-i{x) 
2xHk{x) = Hk+i{x) + 2kHk^i{x), 
which can be derived directly from the Rodrigues formula, to obtain 

.(1) ' 



^l!>{x) = - [Hk+3{x) + 4{k + 3)Hk+i{x) + Ak{k + 3)Hk-i{x)] 
The normalization factor for the harmonic oscillator eigenfunctions is 



(2x2 + 1) 



(2) 



1 



.1/2 



2''k\^ 



so the corresponding factor for the Carinena eigenfunctions is 

1/2 



(1) 



E 



(2) 
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We have therefore exactly solved the Carinena Hamiltonian, and the eigenvalues and 
normalized eigenfunctions are given by 



c 

k 
C, 



X] 



2fc + 3 
2'=+3(fc + 3)!v^ 



'Hu+z{x) + 4(fc + 3)g,.+i(x) + Ak{k + 3)gfc_i(x)] 
(2x2 + 1) 



/2 



where A; = 0, 1, 2 ... To compare with Carinena et al's results we set n = A; + 3 so that 



E. 



c 



XI 



2n-3 

{n - l){n - 2) 
2"n!^/7^ 



nl/2 



[Hn{x) + 4:nHn-2{x) + 4n(n - 3)if„_4(x)] 



(2x2 + 1) 

where n = 3, 4, 5 . . . Note that there is one final state possible, which is the solution of 
^400 = 0, and this will have energy —3 if it exists. Solving this equation gives 

which is of exactly the same form as the previous 0^(x) with n = 0. Even the form of 
the normalisation constant suggested by placing n = in the previous equation turns 
out to be correct, as we show in Appendix A. It follows that the final solution of the 
Carinena Hamiltonian is 



E. 



c 



X 



2n-3 

(n-l)(n-2) 
2"n!x/7r 



1/2 



[Hn{x) + AnHn-2{x) + An{n - 3)g„_4(x)] ^ 
(2x2 + 1) 



-x2/2 



where n = 0, 3, 4, 5 . . . 



We have therefore reproduced the results of Carinena et al in a very direct and 
systematic manner by using the fact that the Cariiiena potential is a partner potential 
of the harmonic oscillator potential. 



4. Partner potentials of the harmonic oscillator 

The solution of the Carinena Hamiltonian in section 3 was motivated by the chance 
observation of one of the authors (JMF) that it is a partner of the harmonic oscillator. 
Let us now turn this around, and ask which Hamiltonians, Hi, are partners to the 
harmonic oscillator, V2(x) = x2, and can thus be solved by the method of section 3. 
The W{x) needed in the factorisation method would then have to satisfy 

'^ + W{xf = V2{x)-\ 
ax 

where we have included an irrelevant constant A. This is an example of a Riccati 
differential equation [6J, first analysed in 1724. The method of solution is to substitute 
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which leads to the equation 

We see that 0(x) is a solution of the original Schrodinger equation, in our case the 
harmonic oscillator. We seem to have gone round in a complete circle, which is not 
surprising since the central idea of the partner potential method is a mapping between 
second order linear and first order non-linear differential equations. 

We also seem to run into a problem in that every solution (f)k{x) of the original 
Schrodinger equation, other than the ground state, will have at least one zero xq. The 
W{x) generated from this will then have a. {x — xo)~^ singularity, and then Vi{x) will 
have a {x — Xo)~^ singularity. This approach would therefore seem to only generate 
soluble potentials which have a singularity at finite x, which might be considered an 
unphysical feature. 

This problem can be overcome in the special case of the harmonic oscillator, for 
which the Schrodinger equation is 

If we set X — iy, this becomes 

which is the original equation with the irrelevant change in constant A — > —A. In other 
words, if we set x — > in the harmonic oscillator eigenf unctions, we get perfectly good 
(t){x) which can then generate W{x) and finally Hi and H2. These (j){x) would not be 
good eigenf unctions as they are not normalizable - they behave like e^^^^ at large x - 
but this is not relevant here. The solutions (j){x) are thus 

(t>p{x) = np{x)e^"'\ 

where the pseudo-Hermite polynomials are given by 

.2 rff 



n,{x) = {-iYH,{ix) = e-^^ — 



e 



They are basically the Hermite polynomials, where the signs of all coefficients are made 
positive. The even solutions 02m (2^) have no real zeros; the odd solutions 02m+i(2^) have 
their only real zero at a; = since they are odd. This leads to a l/x"^ singularity at 
x = 0, which can be regarded as a centrifugal barrier, as in the case of the isotonic 
oscillator. 

Let us examine the first few partner potentials generated in this manner. For A; = 0, 

Wix) = In 00 (2^) = X 
ax 

Vx{x) ^x^-l 

V2{X) = + 1, 
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and this generates the standard ladder operator solution of the harmonic oscillator, since 
this shows that the harmonic oscillator is a partner to itself shifted by 2 units of energy. 
For A; = 1, 



xe 



(f)l{x) 

W(x) = -^\n(f)i(x) = x + - 

ax X 



x'^ + ^ + 1 
x-^ 



Vi{x) 
V2{x) = + 3, 

and the partner potential is an example of the isotonic oscillator. The eigenf unctions of 
the isotonic oscillator generated by this method, 



d 1 
-— +a;+ - 



H^{x)e-^''\ 



are only non-singular for odd m, and hence the energy spectrum is — 2k + 2 where 
/c = l,3,5..., so that the level spacing is double that of the harmonic oscillator. 
For A; = 2, 

ct>2{x) = (4x2 + 2)e^'/2 

d Ax 

w(x) = -i„^,M = x+p^ 

and this is the Carinena potential we have just solved. This is the first new soluble 
potential generated by this method, since the solutions of the harmonic and isotonic 
oscillator are well-known. 
For A; = 3, 

(t>^{x) = (8x=^ + 12x)e^'/2 

4x 



W{x) 



Vi{x) 



x^ 



, i^^' - 3) 
'(2x2 + 3)2 



+ 5 



V2{x) ^x' + 7, 

which is a variant of the Carinena potential which includes a centrifugal barrier term. 
For A; 4, 



24x) 



Mx) = (16x^ + 24x^ + 12)e"'/2 

W(x) = — In 03(0;) ^x + —^, ^ ' , 

^ ^ dx ' (4x4 + 12x2 + 3) 



V,{x) 
V2{x) 



X +16 



(4x4 + 12x2 + 3)^ 



+ 7 



X 



9, 
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which is the harmonic potential plus a rational potential which is regular at the origin, 
and falls off at infinity. 

The common feature of these potentials is that they consist of a harmonic term 
plus an additional rational function which falls off at infinity like a constant times 1/x'^. 
For even p this additional term is symmetric and finite at the origin, and leads to a 
potential which looks hke a harmonic well with an attractive dimple. For odd p the 
additional term is similar to that for even p — 1, plus a centrifugal barrier term. The 
latter means that we consider these potentials in the interval (0, oo), as in the case of 
the isotonic oscillator. 



5. Exact solution of Generalized Carinena potentials 



In this section we will write down the exact closed form solutions of the Generalized 
Carinena potentials, Vq{x), which we define to be the partner potentials of the harmonic 
oscillator generated by the functions (f)p{x) introduced in the last section. The original 
Carinena potential corresponds to the case p — 2. 
We first deduce an expression for Vq{x). Since 



W{x) 

W'{x) 

it follows that 

V2(x) = W^(x) + W'{x) 



d 



dx 

(t)p{x) 



ln0p(x) 



(t)'{x) 



X + 



n'Jx) 



Hp{x) 



(j)p{x) 

(j>p{x)\ Hp{x) 



KM 

Hp{x) 



(t)p{x) 



x^ + 2p+l 



Vi(x) = W^(x) - W'(x) = V2(x) - 2W'(x) = + 2 



n'p{xf-npix)n;{x) 



Hpixy 



+ 2p-l. 



We therefore deduce the formula for the Generalized Carinena potential. 



T^ixf -np{x)n;{x) 



Hpi^xf 



Since 'Hp{x) is a polynomial of order p, we see that the second term in V^(a;) is a rational 
function with numerator a polynomial of degree 2p — 2, and denominator a polynomial 
of degree 2p. As x ^ oo. 



V^{x] 



X'' + —. 

X'^ 



Starting from the original harmonic oscillator eigenvalues, — 2k -\- 1, then 
e'^^ — E^^ — 2(k + p + 1) and E^ — 2k + 3. The unnormahzed eigenf unctions are 



ilj^^\x) = A^ijj^^'^ix) 



d ^ ^ ^ n'pjx) 
dx 'Hp{x) 



2xH,{x) - H'^{x) + ^^^H,{x) 



Hk{x)e-^"'^ 
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n'jx) 



-HJx) 



[np{x)Hk+i{x) + n'^{x)Hk{x) 



-a;2/2 



Hp{x) 



{k + p-i)\ 



The normalization factor is given by 



(2) 



1 



1/2 



If we define n = fc+p+lwe see that the eigenfunctions are given by 



Hpix) 



where 



n 



p + l,p + 2,p + 3. 



the polynomial factors Pn{x) arc related to the Hermite polynomials by 



pi 



[n 



-n- 



■Hn-2i{x), 



and the normalization constant is 



(n — 1) (n — 2) . . . (n — p) 



The corresponding eigenvalues are given by 



1/2 



= -2p + 1 + 2n 



where 



n =p + l,p + 2,p + 3 . . . 



There is one final state possible, which is given by AiIIq^ = 0, and has energy Eq'^ 
so that Eq — —2p + l. This state is then given by 

1 P-^V2 



0, 



(f)p(x) Hpix) 

When p — 2m is even, the normalization constant can be deduced from Appendix A to 
be 

(2m)!2"^^^/^ 



which is exactly what we would obtain if we naively set n = in our equation for the 
normalization constant A^„. It follows that for even p, all the previous formulas are still 
correct, but now we have n = 0,p + l,p + 2,p + 3, . . ., so that we have an equidistant 
level structure but with p levels missing. 

In the case where p is odd, every other level starting with the ground state will have 
a, 1/x singularity at a; = and hence will not be square normalizable, and should be 
eliminated. This gives an equidistant spectrum with twice the spacing of the harmonic 
oscillator levels, exactly as in the case of the isotonic oscillator. 
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6. Summary and Conclusions 

We have obtained the exact solution of Carifiena's quantum non-linear oscillator in an 
economical fashion using the methods of supersymmetric quantum mechanics to show 
that it is a partner of the harmonic oscillator. We have then used this approach to define 
a countably infinite set of generalized Cariiiena potentials V^{x), where p is a positive 
integer, which we then exactly solve in a similar fashion. The case where p is even is 
the most interesting, since all the generated eigenfunctions are normalizable. When p is 
odd, half of the generated eigenfunctions must be removed as they are not normalizable. 

The key technical observation is that the harmonic oscillator eigenvalue equation is 
unchanged under the transformation x ix, so the harmonic oscillator eigenfunctions 
with X replaced by ix can be used to generate the ladder operators A and needed 
to move between supersymmetric partners. More generally one can use any solution of 
a Schrodinger equation to generate the ladder operators; it does not have to be a good 
eigenfunction as normalization is not required. 

We finally note that the isotonic oscillator eigenvalue equation is also unchanged 
under the transformation x ^ ix, and this will be the subject of future investigations. 
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Appendix A. Normalization of ground state eigenfunctions 

In this appendix we evaluate the normalization integral needed for the ground state 
eigenfunctions 

To evaluate this integral define Jp{x) as the indefinite integral 

_a.2 

Jp{x) = J dx = J Gp{xfe~''^dx, 

where for simplicity we have defined Gp{x) — l/7ip{x). We can derive identities for the 
Gp{x) by taking the identities for 7ip{x), 

n'pix) = Hp+iix) - 2xnp(x) 

n'p{x) = 2pnp-i{x), 

and substituting Tlp{x) — 1/Gp{x) to obtain 

Gp{xf = 2xGp{x)Gp+^{x) - G'p{x)Gp+^{x) 

G'p{x)Gp.i{x) = -2pGp{xy. 
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We can now substitute the first identity into the formula for Jp{x) and integrate by 
parts to obtain 



= -Gp{x)Gp+i{x)e~''' + j ([Gp(x)G'p+i(x)]'e~^^ - Gp(x)G'p+i(x)e"^'^) dx 

= —Gp{x)Gp+i{x)e~''^ + j Gp{x)G'p^^{x)e~'^^ dx 

= -G'p(x)G'p+i(x)e~^^ - 2(p+ 1) J Gp+i{xYe~''^ dx 

= -Gp{x)Gp+i{x)e~''^ - 2(p+ l)Jp+i(x), 

where in the fourth step we used the second identity with p replaced by p+ 1. Repeating 
the process and using the formula to evaluate the definite integral gives 

oo 2 oo OO 

Jp{x)\^ = [-Gp{x) + 2{p+l)Gp+2{x)]Gp+,{x)e-^ \^ + 4(p + + 2) Jp+2(x)[ . 

Suppose now that p = 2m is even. The first term on the right hand side clearly 
vanishes as x — ^ oo, but as x ^ the situation is more complicated since G2m+i{x) 
has an 0(1/ x) singularity at x = 0. From the identities for 'Hp{x) we see that 
'^^2m(0) = 2(2m — l)7i2m-2(0) from which we may deduce that 7i2m(0) = (2m)\/m\ 
and hence G2m(0) = m\/{2m)\. The coefficient multiplying G2m+i{x) as x — * is thus 

m\ , , (m + 1)! 

+ 2(2m + 1) ' = 0. 



(2m)! ' '(2m + 2)! 

The 0(l/x) term thus has coefficient zero, and since the next term is 0(x), we see that 
the first term on the right hand side also vanishes as x — 0. If we now multiply by 2 to 
make the region of integration from — oo to oo, we get the result 

4(2m + l)(2m + 2)/2„+2. 
Since Iq = y/n, the result (lA.ll) follows. 

Appendix B. Proof of two Hermite polynomial identities 

In this appendix we prove two identities involving Hermite and pseudo-Hermite 
polynomials, which are need to give an explicit form for eigenfunctions of generalized 
Carinena potentials. 

The first identity we prove is 

H,(.m.) = |:2-^|^^//.,,_..(.). (B,i) 

of which the first few examples are 

1 ■ Hkix) = Hk{x) 

2x ■ Hkix) = Hk+iix) + 2kHk-i{x) 

(4x2 + 2) ■ Hk{x) = Hk+2 + 4(A; + l)Hu{x) + Ak{k - l)Hk-2{x). (B.2) 
To prove this identity we will need the identities 

2xHk{x) = Hk+i{x) + 2kHk-i{x) 
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Hp+i{x) = 2xHp{x) + 2pHp-i{x) , (B.3) 

which are easily derived directly from the relevant Rodrigues formulae. In fact the 
identity fiRTD was originally obtained by using the first of (]B.3|) to derive the examples 
(1B.2I) . and thence spot the general pattern. 

The proof is by induction on p. The case p = is trivially true; the case p = 1 is 
simply the first of flB.3|) . It follows that the formula is true for p = and p = 1. Now 
assume it is true for all values up to p > 2. The formula for p+1 can then be written 

Hp+i{x)Hk{x) = 2xHp{x)Hk{x) + 2pHp-i{x)Hk{x) 

= np{x)Hk+iix) + 2knpix)Hk^i{x) + 2pnp^iix)Hkix), (B.4) 

using the identities ( IB. 31) . The three terms in (IB. 41) take the form 

2pnp^,{x)Hk{x) = 2pE2- -^l ~ 'X Hk+v-i-2:{x) 

^ i\[p - I - i)\ [k + p - I -2i)\ 

t ^'^'-u^iP - ^}r^^^^^^H,^P^^~2^ix), (B.5) 

^ ii{p-iy. {k + p-i-2iy. 

where we were able to extend the summation limit in the first and last terms because 
the additional terms are exactly zero in both cases. Combining the last two terms in 
fIB.SP then yields 

V^oi+i P- ru , -^ik + p-l-i)] 

w^-^ ^""-' hk+p-i- 2.)! ""^-^^-^-^-^"^ 

P'- (fc+P+l-^)! 

" h {i-mp+i~i)\{k+p+i-2i)r'^'^'^''^''^ 

where in the second step we shifted i ^ i — 1, and in the third step we were again able 
to extend the summation limit because the additional term was zero. Adding (]B.6P to 
the first term of (IB.SP gives 

np,,ix)Hkix) -Ej-^^^;^{p+^}f^j^^ 

"yy {k + p + l-^)\ 
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which is just (IB.ip with p replaced by p + 1. This completes the proof by induction. 
The second identity we prove is 

To prove this we note from ( IB. 11) that 
whilst 

n'pix)Hk{x) = 2pnp-i{x)Hk{x) 

- ^\{p -l-^)\{k + p-l- 2z)!^^^-^-^^(^) 

p\ {k + p-l-i)\ 



i=0 



where in the fourth step we shifted z ^ z — 1 , and in the final step we were again able to 
extend the limit of summation. Adding ( IB. 81) and ( IB.OP together then gives the required 
identity (IRT]) . 



References 



[1] Carinena J F, Perelomov A M, Ranada M F and Santander M 2008 J. Phys. A: Math. Gen. 41 
085301. 

[2] Note that we define our Hamiltonian to be twice that used by Carinena, so that our energy values 
are twice his. We use this definition so that the coefficient of the second derivative in Schrodinger's 
equation is -1, which is the standard choice in the factorization and supersymmetry hterature. 

[3] Cooper F, Khare A and Sukhatme U 2001 Supersymmetry in Quantum Mechanics (Singapore: 
World Scientific). 

[4] Infeld L and HuU T E 1951 Rev. Mod. Phys. 23 21-68. 

[5] Fliigge S Practical Quantum Mechanics (Berlin: Springer) p 68. 

[6] Ince E L 1956 Ordinary Differential Equations (New York: Dover) p 23. 



